Impurity scattering in strongly correlated metals close to the Mott transition 



O 

o 



Pascal Lederer^'^ and Marcelo J. Rozenberg^-'^ 

^ Laboratoire de Physique des Solides, CNRS et Universite Paris-Sud, Bdtiment 510, F-91405 Orsay, France. 
^ Departamento de Fsica, UFPE, 50670-901 Recife, Brasil. 
Departamento de Fisica, FCEN, Universidad de Buenos Aires, 
Ciudad Universitaria Pah. I, (1428) Buenos Aires, Argentina. 
(Dated: February 1, 2008) 

This work explores a simple aproximation to describe isolated impurity scattering in a strongly 
correlated metal. The approximation combines conventional one electron scattering theory and 
the Dynamic Mean Field Theory to describe strong correlations in the host. It becomes exact in 
several limits, including those of very weak and very strong impurity potentials. Original electronic 
structure appears at the impurity site when the impurity potential strength is moderate and the 
host is close to the Mott transition. Our results may provide useful guidance for interpretation of 
scanning tunneling microscopy experiments in strongly correlated systems. 
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Introduction. The physics of dilute impurities in broad 
band metals such as Cu or Al is one of the success sto- 
ries of quantum mechanics in the fifties and sixties 
By dilute we mean that physical effect due to interac- 
tions between impurities can be negelected, and that the 
properties described vary linearly with the concentration 
of impurities. For broad band metals, electron-electron 
interactions may be treated as perturbations, because the 
ratio Eel-el /Ek is small compared to unity. Eei-ei ~1- 
2eV is the screened Coulomb interaction between elec- 
trons, and Ek ~10eV is the electronic kinetic energy, 
of the order oi tp, the Fermi energy measured from the 
bottom of the conduction band. Historically, the treat- 
ment of electron-electron interactions in impure metals 
has concentrated for a long time on the Kondo prob- 
lem, for which electron-electron interactions in the host 
where irrelevant More recently, renormalizations of 
the Kondo screening due to correlation effects in the host 
where also considered 0, On the other hand, impu- 
rity scattering with electron-electron interactions on the 
impurity site and in the host, have been treated within 
weak coupling methods in a few simple cases, such as the 
PdNi system [Hj. In recent years, improvements in lo- 
cal probe techniques such as nuclear magnetic resonance 
and scanning tunneling microscopy shed new light on the 
problem of impurity effects in strongly correlated metallic 
materials, in particular with connection to High Tc su- 
perconductivity. For instance, in underdoped supercon- 
ducting cuprates, non magnetic impurities such as Zn or 
Li trigger the appearance of magnetic moments on neigh- 
bouring sites [O, izl . This is a natural phenomenon if the 
two dimensional CuO plane at metallic doping level is 
reasonably well described by some variant of the RVB 
theory @. 

The present paper is dedicated to introduce and discuss 
the results of a simple approach to the problem of iso- 
lated impurity scattering in a strongly correlated metallic 
system close to the Mott transition. 



The understanding of the Mott metal-insulator tran- 
sition in a half filled band has evolved, since the early 
papers by Hubbard with Gutzwiller's treatment of 
strong correlations in narrow band metal [l^, and the 
application of that theory to the half filled band by 
Brinkman and Rice TT]. The latter paper showed that 
close to the critical parameter Uc governing the metal- 
insulator transition, both the Pauli susceptibility x ^nd 
the effective electronic mass m* diverge with a constant 
ratio. This behaviour became understandable thanks 
to the Dynamic Mean Field Theory (DMFT) In 
the simple Hubbard picture of the metal-insulator tran- 
sition, the upper band, with doubly occcupied sites, sep- 
arates at the transition from the lower singly-occupied 
band. Following the DMFT picture, a central peak ap- 
pears between two well separated "Hubbard" bands in 
the metallic phase. The integrated density Z of the cen- 
tral peak goes to zero with the peak width at the tran- 
sition. The peak is analogous to the Kondo resonance 
for the single magnetic impurity problem, except that it 
is a homogeneous property of the system. In addition, 
this narrow feature at the Fermi energy sets a scale for 
a characteristic temperature T*, which goes to zero at 
the metal- insulator transition. Approaching the transi- 
tion, the electronic specific heat coefiicient and the Pauli 
susceptibility diverge at zero T as 1 /T* . 

A number of important properties of the impurity 
in the broad band case are given by the impurity site 
Green's function Goo{t^) = k' ^(^'^' '^)- 
tronic properties of the dilute metal are contained in 
G{k,k',Lu), which is known exactly for the case when 
electron-electron interactions may be neglected alto- 
gether. Goo(w) is entirely computed in terms of the scat- 
tering potential V, which we will assume for simplicity 
to be a point scattering one, and the site diagonal host 
Green's function Gge,{^)- On the other hand, if one con- 
siders the case of a narrow band system, the solution of 
the DMFT equations provides the site diagonal Green's 
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function G'l'^ of a strongly interacting host, where U 
denotes the electron-electron intra atomic "Hubbard" in- 
teraction potential. This Green's function depends solely 
on the momentum independent self energy S(t/, w) and 
on the non- interacting band structure [12| . While the 
DMFT is exact in the limit of large spatial dimensions 
(or lattice connectivity), it can be considered as an ap- 
proximation for a finite dimensional lattice, which is the 
view that we adopt here to describe the host. 

This paper explores the following simple idea: in a cer- 
tain limit, the impurity properties when the host system 
is close to the metal-insulator transition may be approx- 
imated by the substitution of the DMFT host site diag- 
onal Green's function in the expression for the impurity 
site Green's function. This may be considered as a (dy- 
namical) mean field description of the impurity electronic 
properties. This simple approach has several correct lim- 
its, including the case of F — > (and any U) and the 
case [/ ^ (and any V), so one may expect that it may 
be reasonably accurate for the general case. 

The model. We assume that the host is close to the 
Mott transition and has a number of equivalent degener- 
ate bands, such as can be found in a cubic crystal. We 
neglect interband electron-electron interactions. 

The Hamiltonian for each band is then: 

H = Ho + Himp (1) 

Ho = '^Ujcl^^Cj^a + h.c. + U'^ni^^ni^l (2) 



described by 
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Hirap = V{no,] + Uo^l) 



(3) 



The site diagonal pure host Green's function is given 
in the DMFT approximation as ^] : 

G°;,^(^) = ^G°^^(fc,c.) = 5: 1— — (4) 

where e/c is the host electronic dispersion relation in the 
limit [/ = 

From standard isolated impurity scattering theory [l3| , 
the impurity Green's function for [/=0 is: 

The site diagonal impurity Green's function for U =0 
is given by Cr^o = E,,, G{k,k'): 



1-VGM 



(6) 



Scattering potential model. In the scheme that we 
propose here, host correlations are taken into account 
by simply replacing the uncorrelated Green's function 
in Eqini by the DMFT expression of the correlated one. 
The local Green's function at an impurity site is thus 



G^o'^ico) 



(7) 



It is important to realize that this expression has the 
correct limiting behaviours that we mentioned before. In 
fact, one easily checks that in Eql7]the Green's function 
G]^^^ becomes the homogeneous system G'^'^ for 1^=0, 
and the impurity Green's function in a normal (ie, non- 
correlated) host when U=Q. 

It is also worth pointing out that the model ([T][3|) is 
not the exact site diagonal impurity Green's function in 
DMFT scheme In fact, the derivation of the impurity 
Green's function in a correlated host fully within DMFT 
(ie, in the limit of large dimensions), can be obtained in a 
straightforward manner using the "cavity construction" 
|.12i] . However, in that case one realizes that the cavity (ie, 
the Weiss field) of the impurity site exactly coincides with 
the cavity of the clean homogeneous system. Therefore, 
there is no renormalization of the impurity environment 
due to the presence of the scattering potential V, or, in 
other words, the Friedel oscillations on neighboring sites 
are suppressed in the limit of infinite dimensions. Thus, 
that approach is unlikely to provide a correct physical 
description of finite dimensional systems. 

Within our scheme, the poles of G^^ {uj) contain all the 
information about the perturbed local density, the local 
change of the density of states (DOS), and the possible 
occurrence of bound states. In particular, a bound state 
at energy is determined by the simultaneous equa- 
tions: 



l-yRe[G°^^(c^fc)] =0 



Im[G°o^(c.fc)]=0 



(8) 
(9) 



When the first equation is satisfied for cjf,, and 
Im[GQQ^(a;{,)] is small compared to its value in the bulk 
of the pure metal bands, a resonant state is formed, with 

width 



In principle y is a free parameter. However when de- 
scribing actual impurities, the effective scattering poten- 
tial results from the perfect screening constraint. In the 
language of scattering theory, this constraint results in 
the well known Friedel sum rule which connects the 
spherical harmonics phase shifts of the scattered wave 
function to the charge which has to be screened locally 
in order to ensure electrostatic equilibrium of the host 
metal. In the present formulation, the effective scatter- 
ing potential is such that the integrated displaced density 
below the Fermi level must counterbalance the difference 
between the impurity nuclear potential and the host one. 
If the series of 3d transition elements Ti, Cr, Mn, Fe, Co, 
Ni are dissolved as dilute impurities in, say, V2O3, each 
impurity will be described by a potential which attracts 
below or above the Fermi level the number of states re- 
quired for electric neutrality of the alloy. 



CO 

FIG. 1: Top panel: Evolution of the impurity density of states 
p{uj) for increasing V in a non-correlated host (ie, U — 0). 
Bottom panel: Same quantity when the host is in the Mott 
insulator state with f/ = 4. 



Within the present approach, the electronic structure 
around the impurity site in real space can be simply ob- 
tained by Fourier transformation of an expression sim- 
ilar to Eqini where the Fourier transform of ^{i — j), 
G'~'{k), is replaced by the interacting host Green's func- 
tion G'°'^(fc,tj) (cf Eqg]). 

Results. It is now straightforward, using the known 
expression for the host correlated site diagonal Green's 
function provided by the DMFT approach, to find the im- 
purity electronic structure in the whole range of V, U val- 
ues. We use the Iterated Perturbation Theory to solve for 
the host Green's function adopting a semi-circular 
non-interacting DOS. We begin our discussion of our re- 
sults with the systematic study of the local DOS p{lu) 
as a function of the scattering potential V in two ex- 
treme cases. In the upper panel of Fig. [T]we show p{uj) 
for the case of a metallic non-interacting host (standard 
Wolf model). In contrast, in the lower panel we show 
the results for the case of an insulating host which is 
well into the Mott state due to strong correlation effects. 
The local DOS of the Mott insulator corresponds to the 
V = case (continuous line) , where symmetric lower and 
upper Hubbard bands split by a gap « f7 — 2D can be 
observed. The former case describes the well known re- 
sults for a non-interacting metallic host, where the effect 
of the potential is to initially deform the local DOS of 
the host, shifting spectral weight towards lower energies 
(since V < 0), until it becomes strong enough to create a 
bound state out of the conducting band. The condition 
for the formation of the bound state is |F| > 13 for a 
reasonably regular DOS, where D is the half-bandwidth 
of the metal host and is taken equal to one, to set the 
units of the problem. For a Mott insulating host, as V 
increases, spectral weight is also shifted to lower energies 
in both, the upper and the lower Hubbard bands. In ad- 
dition, there is also a transfer of spectral weight from the 
upper Hubard band to the lower one. Interestingly, as the 
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FIG. 2: The evolution of the impurity site p{uj) when the host 
is close to the Mott transition with U = 2.9. The panel series 
show the effect of increasing the strength of the scattering 
potential V. The bottom panel inset shows the detail of the 
evolution of the mid gap resonance for V = 0, 0.5 and 2. 

strength of the scattering potential continues to increase, 
an impurity resonant state forms simultaneously at the 
bottom of the upper and the lower Hubbard bands. The 
threshold value for the appearence of the bound states 
is not as clear cut as in the non-interacting metalic case, 
but remains of the same order of magnitude. It is also 
worth noting that the width of the bound state is very 
narrow for the upper Hubbard band, while substantially 
broader (and more intense) for the lower Hubbard band. 

A correlated metal host is well described with DMFT, 
and the clean (ie V = 0) case is shown in Fig. [5] (top left 
panel). The local DOS shows the characteristic "Kondo" 
peak at uj ~ that corresponds to a renormalized band of 
heavy mass quasiparticles. The effective bandwidth ZD 
of this feature is proportional to its spectral intensity 
(ie the quasiparticle residue Z) and scales with Uc — U, 
where Uc ~ 3-D is the critical value of the interaction 
where a metal-insulator occurs. Spectral weight 1 — .Z is 
thus shifted from low energies to frequencies of order U, 
to form the broad lower and upper Hubbard bands. At 
the transition, the quasiparticle peak becomes infinitely 
narrow and then disapears, leaving behind just the higher 
energy Hubbard bands (cf FigH] lower panel). It is inter- 
esting to observe the systematic effect of the strength of 
the scattering potential on the local DOS for this strongly 
correlated metallic host. As V is stepped up in abso- 
lute value, the DOS is shifted to lower energies in both 
Hubbard bands, as well as in the central peak. At the 
same time a finite amount of spectral weight is trans- 
ferred from the upper Hubbard band to the lower one, 
below the Fermi level, allowing a local screening at the 
impurity site. For \V\ larger than about 1.5-D, a very 
sharp resonance appears at the impurity site. This sharp 
resonance carries a spectral intensity ^ Z, ie the value 
of the quasiparticle residue, and thus can be interpreted 
as a bound state pulled out of the narrow renormalized 
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FIG. 3: Evolution of the impurity density of states p{uj), for a 
fixed strong attractive impurity potential V = 1.5. The panel 
series show the effect of increasing the strength of correlations 
U, as the host is driven across the Mott transition. 



band. This interpretation is further substantiated by the 
results shown in the inset of the lower right panel, that 
shows a detail of the evolution of the DOS near the Fermi 
energy. We observe that the evolution of p{uj) is quahta- 
tively similar as that of the non-interacting metal (Fig[T] 
top panel), but with the frequency scale renormalized 
down by a factor of Z. Eventually, when the strength 
of V is large enough, the quasiparticle band picture near 
the Fermi energy breaks down and a sole resonant state 
at a frequency uj Ri V + U/2 is left. 

In Fig[3] we now consider the case of a strongly at- 
tractive impurity potential \V\ > D, and study the sys- 
tematic evolution of the local DOS as a function of the 
electron-electron interaction parameter [/, as the host 
evolves through a metal-insualtor transition. For weak 
U, the main effect of el-el interactions is to transform 
the impurity bound state to a resonant state, due to the 
tails of the single correlated band. As U is stepped up, 
the two host Hubbard bands form, and the impurity site 
develops a multiple-peak structure. There are two res- 
onances associated with the upper and lower Hubbard 
bands (with the latter being more prominent as V is as- 
sumed attractive), and a central one that is associated to 
the central quasiparticle narrow band. This central fea- 
ture occurs close to the Fermi energy and becomes a very 
sharp resonance at some particular value of C/, and then 
broadens again before disapearing as the host becomes 
a Mott insulator. We find that the spectral strength of 
the mid gap resonance at the impurity site scales roughly 
like the host mid gap quasiparticle peak, and both mid 
gap features vanish together at the Mott transition 

Finally, Fig|4]exhibits the increase of electronic density 
at the impurity site, as a function of \V\ and U. As U 
increases, it requires a stronger attractive potential V to 
screen the same amount of impurity nuclear charge. This 
simply reflects the fact that a large amount of states in 




FIG. 4: The variation of the electronic density at the impurity 
site as a function of for various U values (left), and as a 
function of U for various values of V (right). 

Hilbert space are pushed up in energy away from the 
Fermi level as U increases, so the amount of states avail- 
able to screen he impurity nuclear potential decreases. 

Conclusion. We believe that the simple approximation 
described in this work should be confronted with some 
experimental testing, so as to check its relevance to ac- 
tual strongly correlated impure metals, close to the Mott 
transition. A candidate system to look at would be the 
V2O3 system, with a small concentration c (« 10~^) of 
Sc, Ti, Cr, Mn, Fe, Co, Ni, or Cu. The measurements 
of interest are the transport ones, such as residual re- 
sistivity, temperature dependence of the resistivity, ther- 
moelectric power, etc. and the NMR properties, such as 
Knight shift, Ti. However, perhaps the most clear cut 
validation may be obtained through the direct observa- 
tion of the local electronic density by Scanning Tunneling 
Spectroscopy or precision photoemission measurements. 
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